UNCLASSIFIED 


_ AD  NUMBER _ 

AD917002 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  NOV  1973.  Other 
requests  shall  be  referred  to  Air  Force 
Armament  Laboratory,  Eglin  AFB,  FL  32542. 


_ AUTHORITY 

AFATL  ltr,  14  Nov  1975 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200.20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE. 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


AFATL-TR-73-221 


A  DIGITAL  COMPUTER  PROGRAM  FOR 
EXTRACTING  AERODYNAMIC  COEFFICIENTS 
FROM  SIX-DEGREE-OF-FREEDOM 
DYNAMIC  DATA 


UNIVERSITY  OF  FLORIDA 


TECHNICAL  REPORT  AFATL-TR-7  3-221 
NOVEMBER  1973 


Distribution  limited  to  U.  S.  Government  agencies  only; 
this  report  documents  test  and  evaluation;  distribution 
limitation  applied  November  1973-  Other  requests  for 
this  document  must  be  referred  to  the  Air  Force  A. nament 
Laboratory  (DLMA) ,  ErJ  in  Air  Force  Base,  Florida  32542. 


AIR  FORCE  ARMAMENT  LABORATORY 

All  FOICE  SYSTEMS  COMMAND  •  UNITED  STATES  All  FOICE  I 

>  c 

inn  nfPT 


p. _ 


EGLIN  AIR  FORCE  RASE,  FLORIDA 


A  Digital  Computer  Program  For  Extracting  Aerodynamic 
Coefficients  From  Six-Degree-Of-Freedom  Dynamic  Data 


F.  W.  Stoinbauor 
M.  H.  Clarkson 
T.  E.  Bullock 


Distribution  limited  to  U.  S.  Government  agencies  only; 
this  report  documents  test  and  evaluation;  distribution 
limitation  applied  November  1973.  Other  requests  for 
this  document  must  be  referred  to  the  Air  Force  Armament 
Laboratory  (DLMA),  Eglin  Air  Force  Base,  Florida  32542. 


,'Mv1 


FOREWORD 


This  analysis  was  conducted  by  the  University  of  Florida, 
Cainesville,  Florida,  under  Contract  F08635-73-C-0009,  with  the  Air 
Force  Armament  Laboratory,  Eglin  Air  Force  Base,  Florida.  The  effort 
was  conducted  during  the  period  March  1972  to  March  1973.  Dr.  George 
B.  Findley  (DLMA)  was  program  manager  for  the  Armament  Laboratory. 

This  work  was  partially  supported  by  the  Air  Force  Office  of  Scientific 
Research  (AFOSR)  under  its  project  9871. 

The  principal  investigators  for  the  University  of  Florida  were 
Drs.  T.  E.  Bullock  and  M.  H.  Clarkson. 

This  technical  report  has  been  reviewed  and  is  approved. 


RICHARD  M.  KELLER,1  Colonel,  USAF 
Chief,  Air-to-Surface  Modular 
Guided  Weapons  Division 


ABSTRACT 


The  development  of  a  digital  computer  program  to  extract  aero¬ 
dynamic  coefficients  from  dynamic  data  for  six-degree-of-freedom 
systems  is  presented.  The  derivation  of  a  system  mathematical  model 
is  discussed  in  detail.  Results  and  associated  problems  of  extracting 
coefficients  from  one-,  two-,  three-,and  six-degree-of-freedom  systems 
data  are  also  presented. 
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SECTION  I 


INTRODUCTION 


Early  methods  of  extracting  aerodynamic  coefficients  from  dynamic 
data  required  assumptions  and  limitations  so  that  the  equations  of 
motion  could  be  solved  in  closed  form.  Therefore,  usually  only  linear 
aerodynamics  were  considered.  As  a  result,  the  range  of  motions  and 
the  number  of  coefficients  extracted  were  severely  limited. 

The  method  of  extracting  coefficients  by  means  of  parametric  differ¬ 
entiation  developed  by  Chapman  and  KirkXl)  is  not  restricted  by  the 
requirement  of  linear  aerodynamics.  In  this  report  the  method  of  para¬ 
metric  differentiation  is  used  to  develop  a  six-degree-of-freedom  digital 
computer  program  to  extract  aerodynamic  coefficients  from  free  flight 
data.  The  program  is  an  extension  of  the  one-  and  three-degree-of-freedom 
programs  of  Daniel  and  Bullock(2),  respectively,  and  draws  on  their 
experience  in  developing  those  programs. 

The  equations  of  motion  for  the  six-degree-of-freedom  mathematical 
model  were  developed  so  that  the  aerodynamic  coefficients  presented  by 
Holmes (3)  may  be  used.  The  model,  which  is  intended  for  use  in  an  aero- 
ballistic  wind  tunnel  range,  uses  a  fixed  plane  axis  system  to  represent 
the  angular  orientation  of  the  body  with  respect  to  a  tunnel  fixed  axis 
system.  By  definition,  the  fixed  plane  axes  are  free  to  pitch  and  yaw 
with  the  body  but  do  not  roll  with  the  body.  The  relationship  between 
the  axis  systems  is  depicted  in  the  diagram  below. 


),p  Tunnel  Fixed  axes 
)  Fixed  plane  axes 
)B  Body  fixed  axes 
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SECTION  II 
EQUATIONS  OF  MOTION 


The  derivation  of  the  equations  of  motion  for  the  six-degree-of- 
freedom  model  used  in  the  program  assumes  that  the  missile  is  regarded 
as  a  rigid  axisymmetric  body  moving  with  velocity  Vx  relative  to  a  wind 
tunnel  axis  system.  In  addition,  the  body  fixed  axes  are  chosen  to 
coincide  with  the  principal  axes  of  the  missile. 

The  equations  for  translational  and  angular  motion,  based  on  Newton's 
second  law,  may  be  written  as 

"  3F  5t  ’  'ft  <» 

and 

dF  +  “fp  x  ~h  =  ”fp  (2) 

is  the  resultant  external  force. 

is  the  angular  velocity  of  the  fixed  plane  axes  with 
respect  to  the  tunnel  fixed  axes. 

is  the  moment  of  momentum. 

is  the  resultant  external  moment. 


where : 


FP 


^Vp 


The  equations  of  motion  above  provide  a  form  suitable  for  fitting 
to  the  data.  First,  however,  it  is  necessary  to  define  the  orientation 
of  the  fixed  plane  axes  with  respect  to  the  aerodynamic  data  axes,  which 
contain  the  cameras  that  recorded  the  motion  and  position  of  the  body 
during  flight,  and  then  to  define  the  tunnel  axes  (assumed  inertial)  with 
respect  to  the  fixed  plane  axes. 


Choosing  the  body  fixed  axes  to  lie  along  the  principal  axes  of  the 
missile  results  in  the  products  of  inertia  being  zero.  Thus,  the  angular 
momentum  vector  may  be  expressed  in  terms  of  the  angular  velocity  and  the 
moments  of  inertia. 


h=hi  +  hj+hk=Iui  +  Iwj+It)k 
x  y  z  xx  y  y  z  z 


(3) 
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Recalling  the  relationship  between  the  fixed  plane  axes  and  the 
body  fixed  axes,  equation  (2)  may  be  written  in  component  form 
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Recalling  that  the  body  fixed  axes  were  principal  axes  implies  that 

I  *  I  -  I  (5) 

y  z 

Now  expressing  the  angular  velocity  components  of  the  fixed  plane  axes 
in  terms  of  the  Euler  angles  yields 

A 

id  m  -  ip  Sin  0  ■  p 
XFP 

u)  ■  0  =  q  (6) 

yFP 

a)  ■  i|i  Cos  0  ■  r 
ZFP 


(7) 
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and  the  angular  velocity  components  of  the  body  fixed  axes  are 

ux  “  ♦  ~  'P  Sin  0  ■  p 

u)*0  *  q 

y 
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which  have  time  derivatives 
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Now  applying  equations  (5),  (6),  (7)  and  (8)  to  the  first  of  equation  (4) 
results  in 

Xx  L  ^  _  ^  Sin  0  -  ije  Cos  (9) 

FP 

Similarly  the  second  of  equations  (4)  becomes 

10  +  Cos  6  Ix  Sin  e  J  +4 1  Sin  0  1^  Cos  6  ■  My  (10) 

Rearranging  terms  yields 


10  +  Cos  6  l^p  +  I  ij  Sin  6  * 
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or 
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Now  operating  in  the  same  manner  on  the  third  of  equations  (4)  yields 
I  |^i|)  Cos  0  -  ij)0  Sin  0^  -  i}>  Sin  0I0-0Ix^-i|»  Sin  6^  = 


FP 

(13) 


Rearranging  terms  yields 

I  Cos  0  -  i|)0  Sin  0 
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Consider,  now,  equation  (1)  for  translational  motion, 
written,  directly,  in  component  form  as 


(14) 


(15) 


It  may  be 
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The  definitions  of  the  resultant  aerodynamic  forces  above  and  the 
resultant  aerodynamic  moments  were  represented  in  terms  of  resultant 
aerodynamic  force  and  moment  coefficients,  Cy,  Cy,  Cg  and  C^>  C^, 
which  lie  along  the  aerodynamic  data  axes.  For  the  translational  equations 
of  motion  it  was  first  necessary  to  prescribe  how  the  components  of  each 
aerodynamic  coefficient  along  the  tunnel  fixed  axes  would  be  determined 
in  terms  of  the  fixed  plane  axes.  Then  for  all  of  the  equations  of  motion 
it  was  necessary  to  transform  the  components  along  the  fixed  plane  axes 
in  terms  of  the  aerodynamic  axes.  The  transformations  for  the  transla¬ 
tional  equations  of  motion  were 
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The  first  and  second  transformations  are  straightforward  transfor¬ 
mations  through  the  Euler  angles  \l>  and  6  from  the  tunnel  fixed  axes  to 
the  fixed  plane  axes.  The  third  transformation  is  the  transformation 
of  the  coefficients  in  terms  of  the  fixed  plane  axes  to  the  coefficients 
in  terms  of  the  aerodynamic  data  axes.  The  application  of  the  transfor¬ 
mation  matrices  yields  the  equations 


L(*)  1.(0) 


(20) 
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or,  in  expanded  form 


X_  *  \  C  (Cos  0  Cos  ij/)  -  C  (Sin  ^  Cos  $  + 
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Before  arriving  at  the  final  form  of  the  angular  equations  of 
motion  a  transformation  to  obtain  the  components  of  the  coefficients 
along  the  aerodynamic  axes  employing  the  third  transformation,  equation 
(19),  mu6t  be  carried  out.  The  resulting  equations  are 
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or,  in  expanded  form, 
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0  ■  -  (p  +  ip  Sin  0)  ip  Cob  0  +  (C^  Cos  $  + 


C„  Sin  ip  + 
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V$” 
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The  resultant  coefficients  in  the  equations  of  motion  are  defined  in 
the  following  fashion: 

C  -  C  +  C  a2 
*  xo  xa2 
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Individual  coefficients  are  defined  in  the  list  of  symbols. 

In  order  to  avoid  ambiguities  which  might  occur,  the  Euler  angles 
ip,  6  and  p  are  limited  to  the  following  ranges: 
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For  non-planar  motion,  that  is,  pitching  and  yawing  motions  occurring 
simultaneously,  the  limits  on  the  ranges  of  the  Euler  angles  and  0 
should  be 


tt 
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to  obtain  reasonable  accuracy  of  the  coefficients  extracted  without  an 
excessive  number  of  iterations. 
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SECTION  III 


METHOD  OF  EXTRACTING  COEFFICIENTS 
AND  DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


1.  Chapman  and  Kirk  Coefficient  Extraction  Method 

The  value  of  using  parametric  influence  coefficients  in  the  analysis 
of  dynamic  systems  has  been  recognized  for  some  time.  The  following 
briefly  describes  the  general  scheme  developed  by  Chapman  and  Kirk  to  use 
the  method  of  parametric  influence  coefficients  for  determining  aerody¬ 
namic  coefficients.  A  more  detailed  presentation  of  the  theory  is  given 
in  references  1  and  2. 

The  basis  of  the  method  is  the  minimization  of  the  deviations  of  a 
set  of  experimental  data  from  a  calculated  motion.  The  system  model 
that  yields  the  calculated  motion  is  given  by  the  set  of  differential 
equations 


x  =  f (x,C»t)  ,  x(0)  *  a  (29) 

where  x(t)  is  an  (n  x  1)  state  vector,  f  is  the  (n  x  1)  vector-valued 
function,  and  a  is  the  (n  x  1)  vector  of  initial  conditions.  The  set  of 
experimental  data,  z(t),  are  the  components  of  the  state  vector.  Assuming 
that  x(t)  is  measured  for  0  _<t  ^x,  then  it  is  desired  to  find  the  para¬ 
meters  C  which  minimize  the  expression 

T 

MSQE  *  7  /'  (  x(t)  "  z(t)  F  (^(t)  {  ;:(t)  -  z(t)  }  dt  (30) 

o 

where  z(t)  are  the  experimental  data  corresponding  to  the  calculated  motion 
of  the  state  vector  x(t)  and  Q«(t)  is  an  (n  x  n)  weighting  matrix  whose 
purpose  is  to  give  weight,  or  value,  to  only  those  components  of  the  state 
vector  x(t)  for  which  measured  experimental  data  are  available. 

The  method  used  to  determine  the  parameters  C  that  satisfy  equation 
(30)  was  an  Iterative  one.  For  each  iteration  a  calculated  motion  and  a 
corresponding  mean  square  error  (MSQE)  was  determined.  If  the  change  in 
the  root  of  the  mean  square  error  was  not  less  than  a  predetermined  value, 
the  parameters  c  were  updated,  or  corrected,  toward  that  end.  This  was 
accomplished  by  Integrating  the  set  of  equations  obtained  by  taking  the 
partial  derivatives  of  each  of  the  equations  of  motion  with  respect  to 
the  parameters  of  interest.  These  will  be  referred  to  as  parametric 
differential  equations  in  this  paper.  The  solutions  of  the  parametric 
differential  equations,  parameter  influence  coefficients,  were  then  used 
to  construct  the  (p  x  p)  matrix  of  what  will  be  referred  to  as  parametric 
influence  coefficients 
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For  example,  if  there  were  4  initial  conditions  and  6  coefficients, 
or  10  parameters  of  interest,  sixty  second  order  parametric  differ¬ 
ential  equations  were  integrated  to  obtain  the  10  x  10  matrix  of  parametric 
influence  coefficients.  Simultaneously,  the  gradient 


3  f 

3  C , 


3  f 

3  C, 


(31) 


(32) 


was  obtained.  Then  the  p  x  1  matrix  of  parameter  corrections,  AC,  was 
found  by 

[AC]  =  [Af1  [B]  (33) 


The  parameters  for  the  next,  or  £  +  1,  iteration  were 


lCb  +  1  *  1°]*  +  <4C)t 


(34) 


Once  the  parameters  were  corrected,  a  new  calculated  motion  was 
determined  by  integrating  the  equations  of  motion.  The  entire  process 
was  repeated  until  the  predetermined  value  for  the  change  of  the  root 
of  the  mean  square  error  was  satisfied,  and  the  process  was  said  to  have 
converged,  or  the  maximum  number  of  iterations  allowed  was  execeeded  and 
the  program  was  terminated. 

2.  Development  of  the  Program 

The  program  is  written  in  Fortran  IV  for  use  primarily  on  an  IBM 
360/65  or  370/165  computer.  The  program  provides  the  user  with  three 
general  options: 

1.  Flight  simulation 

2.  Coefficient  extraction 

3.  Flight  simulation  with  punched  output  of  state  vector 
component  time  histories. 

The  paragraphs  that  follow  describe  the  functions  of  the  main  program, 
its  subroutines,  and  the  program  options.  A  flow  chart  and  a  complete 
listing  of  the  program  and  the  required  data  input  form  are  given  in  the 
appendices. 


1.0 


The  function  of  the  main  program  is  to  control  the  flow  of  the 
program  in  accordance  with  the  options  chosen.  To  do  this,  t lie  main 
program  reads  and  writes  all  input  and  output  Information,  organizes 
the  information,  and  calls  the  subroutines  to  use  it.  The  main  program 
does  all  of  the  calculations  necessary  to  determine  if  convergence  has 
been  achieved  and  all  of  the  calculations  preparing  for  each  iteration. 

Subroutine  ADDUM  integrates  the  equations  of  motion  and  the  para¬ 
metric  differential  equations.  The  numerical  method  used  is  a  fourth 
order  Runge-Kutta  starter  solution  and  a  fourth  order  Adams-Bashforth 
predictor-corrector  method  for  integrating.  This  subroutine  is  des¬ 
cribed  in  detail  in  reference  5. 

Subroutine  XD0T1  computes  current  values  of  the  derivatives  of 
the  set  of  first  order  equations  to  which  the  equations  of  motion  have 
been  reduced  as  required  by  ADDUM.  The  subroutine  also  computes  the 
value  of  the  derivative  of  the  mean  square  error,  which  is  Integrated 
simultaneously  with  the  equations  of  motion  by  ADDUM  when  the  coefficient 
extraction  option  is  specified. 

Subroutine  0UT1  stores  the  results  of  the  numerically  integrated 
equations  of  motion  during  each  iteration  until  convergence  is  tested. 

Subroutine  XD0T2  computes  current  values  of  the  derivatives  of 
the  set  of  first  order  equations  to  which  the  parametric  differential 
equations  have  been  reduced  as  required  by  ADDUM. 

Subroutine  0UT2  calculates  the  elements  of  the  parametric  influence 
coefficient  matrix,  [A],  and  the  state  vector  difference  matrix,  [B],  for 
use  in  the  main  program. 

Subroutine  MINV  inverts  the  p  x  p  matrix  of  parametric  influence 
coefficients  using  a  standard  Gauss-Jordan  method  and  is  described  in 
detail  in  reference  6. 

Subroutine  PL0T9  is  a  printer-plotter  routine  intended  to  give  the 
program  user  a  visual  understanding  of  the  angular  orientation  of  the 
missile  as  calculated  by  the  equations  of  motion. 

The  amount  of  input  data  required  by  the  program  is  determined 
by  the  program  option  chosen.  The  specific  data  in  each  option  are 
delineated  in  the  following  paragraphs.  The  formats  and  units  of 
entries  on  specific  data  cards  may  be  found  in  Appendix  III. 

(1)  Flight  simulation 

(a)  Program  control  codes 

These  integer  constants  tell  the  program  which  program  options  are 
in  effect  and  which  equations  of  motion  are  to  have  values  computed  fcr 
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their  derivatives  in  XD0T1.  The  integrated  values  of  all  other  equations 
of  motion  are  set  to  zero.  The  purpose  of  allowing  the  program  user  to 
specify  the  equations  of  motion  that  will  have  nonzero-integrated  values 
is  to  avoid  unnecessary  computation,  thus  reducing  execution  time. 

(b)  Integration  constants 

The  integration  constants  include  the  numerical  integration  step 
size,  frequency  of  storage  of  integrated  values,  time  at  which  integration 
is  to  stop,  initial  time,  YES  or  NO  codes  to  specify  printer  plots  of 
each  of  the  three  angular  motions,  and  the  number  of  fins  on  the  missile. 
The  number  of  fins  choice  allows  the  user  to  specify  a  four- finned 
missile  or  an  unfinned  projectile. 

(c)  Aerodynamic  and  physical  constants 

The  aerodynamic  constants  are  air  density  and  the  free  stream 
velocity  that  are  specified  during  the  flight  simulation.  The  physical 
constants  are  the  body  cross-sectional  area  (neglecting  fins),  body 
diameter  or  equivalent,  spin  rate  of  the  body  at  time  zero,  gravitational 
acceleration  due  to  the  earth,  moment  of  inertia  about  the  longitudinal 
axis,  moment  of  inertia  about  the  axes  normal  to  the  pitch  and  yaw  planes, 
and  mass  of  the  body. 

(d)  Aerodynamic  coefficients 

For  flight  simulation  the  aerodynamic  coefficients  values  are  con¬ 
stant  and  are  not  altered  by  the  program. 

(e)  Initial  conditions 

These  values  are  the  initial  conditions  for  the  equations  of  motion. 
Like  the  aerodynamic  coefficients,  they  are  constant  and  are  not  altered 
by  the  program. 


(f)  Printer  plotter  constants 

These  constants  are  required  only  if  the  plot  option  was  specified 
in  the  integration  constants,  The  constants  are  the  width  of  the  plot, 
value  of  the  initial  point,  type  of  plot,  field  type  for  the  data  point 
values  printed,  and  a  scale  factor. 

(2)  Coefficient  extraction 

(a)  Output  labels 

These  labels  allow  the  program  to  identify  the  extracted  values 

and  the  estimated  standard  deviations  with  appropriate  labels. 

(b)  Program  control  codes 


In  addition  to  these  listed  for  flight  simulation,  there  are  con¬ 
stants  to  specify  initial  conditions  and  aerodynamic  coefficients  to 
be  adjusted,  values  for  the  weight  factors  in  equation  (30),  maximum 
number  of  iterations  allowed, and  convergence  tolerance  before  the  iteration 
process  is  automatically  terminated. 

(c)  Data 

The  initial  condition  and  coefficient  values  input  are  now  guesses 
and  not  constant  values.  In  addition! values  for  the  experimental  data 
points  of  the  state  vector  components  must  also  be  read. 

(3)  Flight  simulation  with  punched  output 

The  input  for  this  option  differs  from  the  flight  simulation  only 
in  the  addition  of  a  program  control  code  to  specify  which  state  vector 
components  are  to  be  punched  on  cards. 
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SECTION  IV 


RESULTS  OF  TEST  CASES 


Eleven  test  cases  of  the  program  were  run  to  check  its  operation. 

The  test  cases  began  with  a  one-degree-of- freedom  case  and  were  increased 
to  a  six-degree-of-freedom  case.  In  all  but  two  cases,  the  initial  con¬ 
ditions  and  aerodynamic  coefficients  used  to  generate  the  data  for  the 
extraction  program  were  known.  This  provided  the  easiest  method  for 
checking  the  validity  of  extracted  initial  conditions  and  coefficients. 

In  the  two  cases  where  initial  conditions  and  coefficients  were  not 
known,  the  extracted  values  were  compared  with  those  obtained  from  the 
same  data  by  Daniel  using  UFPLANAR.  The  reason  for  investigating  these 
two  cases  was  to  check  the  capability  of  the  program  to  handle  noisy 
data.  The  noise  was  simulated  by  random  measurement  errors  in  UFN0ISE(2). 
A  table  of  the  results  of  the  eleven  cases  may  be  found  in  Appendix  IV. 

The  two  cases  with  noisy  data  considered  one-degree-of-freedom 
cases  with  linear  and  non-linear  static  restoring  moment  and  pitch 
damping  coefficients.  As  intuition  would  lead  one  to  expect,  the  es¬ 
timated  standard  deviations  of  the  values  extracted  from  noisy  data 
were  much  larger  than  the  standard  deviations  of  the  values  extracted 
from  data  without  noise.  The  standard  deviations  of  the  values  extracted 
from  data  without  noise  were  essentially  zero,  as  they  should  have  been, 
since  the  data  were  generated  from  the  same  equations  of  motion.  How¬ 
ever,  the  important  result  was  that  the  number  of  iterations  required 
for  convergence  was  the  same  for  both  types  of  cases.  This  is  very 
desirable  from  a  computing  standpoint  because  free  flight  test  data  will 
most  certainly  be  noisy. 

As  stated  previously  in  Chapter  II,  the  mathematical  model  of  the 
missile  was  restricted  to  low  angles  of  attack  for  multiangular  degree 
of  freedom  cases.  In  order  to  quantitatively  demonstrate  the  necessity 
for  this  restriction,  two  cases  were  run  with  initial  pitch  and  yaw 
angles  both  equal  to  20  degrees  in  the  first  case  and  30  degrees  in  the 
second.  The  20-degree  case  required  a  reasonable  six  iterations  to 
extract  initial  conditions  and  coefficients.  On  the  other  hand,  the  30- 
degree  case  required  eleven  iterations  to  extract  the  correct  values. 

Several  multi-degree -of-freedora  combinations  of  angular  and  trans¬ 
lational  motions  were  among  the  cases  run.  It  was  found  for  these  cases 
that  the  extraction  process  had  to  be  a  two-  or  three-step  process, 
depending  on  the  complexity  of  the  case.  The  necessity  for  this  proce¬ 
dure  is  a  matter  of  the  relative  sensitivity  of  the  parameters.  This 
sensitivity  may  be  observed  by  comparing  the  magnitudes  of  the  elements 
along  the  main  diagonal  of  the  influence  coefficient  matrix.  If  a 
parameter  is  either  insensitive  or  too  sensitive  to  the  motion  of  the 
missile,  it  will  cause  the  adjustment  of  the  parameters  from  iteration 
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to  iteration  to  be  incorrect,  that  is,  too  small  or  too  large.  By 
carrying  out  the  extraction  process  in  a  certain  order  of  steps,  this 
problem  can  be  avoided.  The  steps  should  be  as  follows: 

1.  Extract  initial  conditions  and  coefficients  related 
solely  to  translational  motion. 

2.  Extract  initial  conditions  and  coefficients  related 
solely  to  angular  motion. 

3.  Extract  coefficients  related  to  interacting  motions,  such 
as  the  magnus  forces  and  moments. 

The  order  of  the  steps  is  as  important  as  the  steps  themselves. 

The  translational  motion  must  be  dealt  with  first  since  the  formulation 
of  the  total  angle  of  attack  requires  the  inclusion  of  the  velocity 
components. 

For  cases  considering  only  pitching  and/or  yawing  motions  with  a 
rolling  motion,  the  extraction  process  requires  only  two  steps: 

1.  Extract  initial  conditions  and  coefficients 
related  solely  to  rolling  motion. 

2.  Extract  initial  conditions  and  coefficients 
related  to  pitching  and/or  yawing  motion. 

It  should  be  noted  that  the  two  aforementioned  processes  are  recom¬ 
mendations. 
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SECTION  V 


CONCLUDING  REMARKS 


In  summary,  the  purpose  of  this  report  was  to  construct  a  six- 
degree-of- freedom  digital  computer  program  which  extracts  aerodynamic 
coefficients  from  free  flight  test  data  using  the  Chapman  and  Kirk 
scheme.  The  mathematical  model  chosen  for  the  program  is  somewhat 
arbitrary;  the  model  has  limitations  such  as  the  number  and  type  of 
aerodynamic  coefficients  and  magnitude  of  the  angles  of  attack  for 
multiangular  degree  of  freedom  cases  as  has  already  been  shown.  These 
limitations  are  not  a  function  of  the  extraction  scheme.  With  this  in 
mind,  the  program  was  designed  to  be  readily  adaptable  to  a  wide  range 
of  mathematical  models.  Major  portions  of  the  model  are  incorporated 
in  subroutines  to  facilitate  any  changes.  For  instance,  the  input  and 
output  of  information  and  the  associated  operations  are  contained  in  the 
main  program,  the  model  equations  of  motion  are  in  XDOTl,  and  the  para¬ 
metric  differential  equations  are  in  XD0T2.  In  addition  to  making 
program  changes  an  easier  process,  this  feature  allows  major  portions 
of  the  program  to  be  bypassed,  depending  on  the  program  option  chosen. 

Segmenting  the  program  into  subroutines  and  using  only  those 
necessary  in  a  given  run  is  a  method  of  keeping  execution  time  to  a 
minimum.  However,  since  the  program  was  designed  to  be  capable  of 
handling  a  maximum  of  six  degrees  of  freedom,  its  operation  on  only 
one  or  two  degrees  of  freedom  is  relatively  costly.  Thus,  for  maximum 
efficiency  its  use  should  be  limited  to  multi-degree  of  freedom  cases. 

There  are  several  areas  of  the  program  to  be  considered  for  further 
study  or  refinement.  The  first  is  the  system  model.  One  can  see  that 
this  is  an  area  of  problem  trade-offs.  A  very  general  model  capable 
of  handling  a  greater  and  more  varied  number  of  aerodynamic  coefficients 
is  more  desirable  from  a  purist's  standpoint.  However,  the  increased 
complexity  and  execution  time  of  such  a  model  is  undesirable.  Of  course, 
the  model  may  be  designed  to  satisfy  only  certain  requirements  and  yield 
good  execution  times  but  at  the  expense  of  generality.  In  addition, 
important  aspects  of  any  change  are  the  time  and  effort  necessary  to 
make  that  change. 

A  second  area  that  should  be  considered  is  the  programming  techniques 
used  in  constructing  the  program.  The  program  was  written  in  a  straight¬ 
forward  manner  rather  like  translating  English  to  a  foreign  language  word 
by  word  to  make  the  program  logic  more  understandable  to  the  user. 

Although  efficient  programming  techniques  would  reduce  execution  times, 
this  is  accompanied  by  a  program  that  would  be  less  understandable  to  the 
user. 


Other  areas  that  might  be  considered  are  much  more  complex.  From 
informal  discussions  with  Chapman  and  other  sources,  such  as  Meissinger , 
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the  author  feels  that  the  parameter  influence  coefficients  and  para¬ 
metric  influence  coefficients  in  the  [A]  matrix  are  another  important 
area  for  further  study.  It  has  already  been  found  that  the  elements 
of  the  [A]  matrix  can  be  an  important  guide  to  the  sensitivity  of  a 
parameter  to  the  motion  of  the  system  model.  With  further  study  it 
might  be  possible  to  determine  not  only  the  numerical  value  of  a  para¬ 
meter,  but  also  its  importance  to  the  system  relative  to  the  other 
parameters  in  a  quantitative  sense  rather  than  just  a  qualitative  one. 
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APPENDIX  I 


PROGRAM  VARIABLE 
H 

ITO 

TMAX 

TZERO 

XZO(I) 

STDIC(I) 

COEF(I) 

STDC(I) 

ICADJ(I) 

CADJ(I) 

QW(I) 

MAXIT 

TOL 

NPTS 

N 

RO 


SIX-DEGREE-OF-FREEDOM  NOMENCLATURE  LIST 
AND  PROGRAM  LISTING 


Nomenclature  list  (partial) 

MATH  SYMBOL  DEFINITION 


At 


t 

max 


t 

o 


v*> 


p 


Numerical  integration  step  size  (sec) 

Frequency  of  numerical  integration 
output 

Cutoff  time  for  numerical  integration 
(sec) 

Initial  time  for  numerical  integration 
(sec) 

Initial  condition  labels 

Initial  condition  standard  deviation 
labels 

Coefficient  labels 

Coefficient  standard  deviation  labels 
Initial  conditions  extracted  code 
Coefficients  extracted  code 
Weight  factor 

Maximum  number  of  iterations  before 
program  terminates 

Convergence  criteria  for  change  in 
root  mean  square  error 

Number  of  experimental  data  points 

Nunber  of  first  order  differential 
equations 

Air  density 
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PROGRAM  VARIABLE 


MATH  SYMBOL 


DEFINITION 


V 

V 

Wind  tunnel  velocity  (ft/sec) 

AR 

A 

Body  reference  area  (ft) 

D 

d 

Body  reference  diameter  (ft) 

P 

P 

Body  spin  rate  (rad/sec) 

G 

g 

Gravitational  constant  (ft/sec) 

AIX 

I 

X 

Moment  of  inertia  about  an  axis 
longitudinally  through  the  CG 
of  the  body  (slug  ft) 

AM 

m 

Mass  of  the  body  (slugs) 

CLA 

C£a 

Rolling  moment  coefficient  (rad) 

CLP 

Clp 

Roll  damping  coefficient  (rad) 

CMA 

C  - 
ma 

Static  pitching  moment  coefficient  (rad) 

CNA 

C  - 
not 

Pitching  moment  (due  to  fins) 
coefficient  (rad) 

CNPA 

C  - 
npa 

Magnus  moment  coefficient  (rad) 

CMQO 

C 

mqo 

Pitch  and/or  yaw  damping  coefficient 
(rad) 

CXO 

C 

xo 

Drag  coefficient  (rad) 

CYA 

C  - 

ya 

Side  force  coefficient  (rad) 

CYAP 

C  - 

yap 

Magnus  force  coefficient  (rad) 

CZA 

c  - 

za 

Normal  force  coefficient  (rad) 

<  >2 

(  )2 

Second  order  term 

t  >3 

<  >3 

Third  order  term 

Angle  of  pitch  in  X-Z  plane  (rad) 
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PROGRAM  VARIABLE 


MATH  SYMBOL 


DEFINITION 


X 

Y 

Z 

(  ) 

(  ) 

DATUM(I) 

DCALC(I) 

IEQ(I) 

IP 

IT 

IFE 

NF 


x 


y 


z 


d_ 

dt 


z(t) 

x(t) 


Angle  of  yaw  in  X-Y  plane  (rad) 

Angle  of  body  roll  relative  to 
fixed  plane  axis  system  (rad) 

X  position  of  body  relative  to 
tunnel  reference  point  (ft) 

Y  position  of  body  relative  to 
tunnel  reference  point  (ft) 

Z  position  of  body  relative  to 
tunnel  reference  point  (ft) 

First  derivative  with  respect 
to  time 

Second  derivative  with  respect 
to  time 

Experimental  data  point  values 

Calculated  point  values 

Equations  of  motion  to  be  integrated 
code 

Psi  plot  code 
Theta  plot  code 
Phi  plot  code 
Number  of  fins  code 
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RETURN 

END 


APPENDIX  II 


FLOW  CHART  OF  COMPUTER  PROGRAM 


Read  input  codes 


Flight  simulation 

I 

Read  and  v;rite  initial 
conditions  and  coefficients 
and  other  data 


- : — 1 - 

Integrate  equations  of  motion 


Coefficient  extraction 

I 

Read  and  write  initial 
conditions  and  coefficients 
and  other  data 

J 

Read  time  history  data 

_ I 


1 - 

Write  or  punch  flight 
simulation  time  history 

J 

End 


\ - 

Convergence  achieved. 

Write  extracted  parameters 

I 

End 


Calculate  mean  square  error 
and  test  root  mean  square 
error  for  convergence 


Convergence  failed.  Maximum 
number  of  iterations  exceeded 


End 


Calculate  parameter  corrections 

I 

Update  parameters 

I 

Go  to  integration  of 
equations  of  motion 
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APPENDIX  III 
DATA  INPUT  FORMAT 


CASE  I  Flight  simulation 


CARD 

COLUMNS 

FIELD 

1 

1 

11 

2 

1-6 

11 

explanation 

Mode  =  0 

IEQ(I)  Equations  of  motion  to  be  inte 
grated.  In  ordei;  the  equations  are: 

Yawing 
Pitching  (6) 

Rolling  (4>) 

X  translation  (X) 

Y  translation  (Y) 

Z  translation  (Z) 


3 


4 


Integration  constants 


1-5 

F5.3 

H 

6-8 

13 

ITO 

9-13 

F5.3 

TMAX 

14-18 

F5.3 

TZERO 

19 

11 

IP  0 :no,  l:yes 

20 

11 

IT  0 : no ,  l:yes 

21 

11 

IFE  0:no,  l:yes 

22 

.11 

NF 

Aerodynamic  constants 

1-11 

F11.8 

RO 

12-22 

F11.4 

V 
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CARD 

COLUMNS 

FIELD 

EXPLANATION 

23-33 

F11.4 

AR 

34-44 

F11.4 

D 

45-55 

F11.4 

P 

5 

Aerodynamic  constants 

1-11 

F11.8 

g 

12-22 

F11.4 

A IX 

13-33 

F11.4 

AI 

34-44 

F11.4 

AM 

6 

Aerodynamic  coefficient  values 

1-10 

F10.4 

CLA 

11-20 

F10.4 

CLA3 

21-30 

F10.4 

CLP 

31-40 

F10.4 

CMA 

41-50 

F10.4 

CMA3 

51-60 

F10.4 

CNA 

61-70 

F10.4 

CNA3 

71-80 

F10.4 

CNPA 

7 

Aerodynamic  coefficient  values 

1-10 

F10.4 

CNPA3 

11-20 

F10.4 

CMQO 

21-30 

F10.4 

CMQ2 

31-40 

F10.4 

CX0 

41-50 

F10.4 

CXA2 

51-60 

F10.4 

CYA 

61-70 

F10.4 

CYA3 
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CARD 


8 


9 


10 


11 

(12,13) 


COLUMNS 

71-80 

1-10 

11-20 

21-30 

1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 
71-80 

1-10 
11-20 
21-30 
31-40 
41-50 


FIELD 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

F10.4 


EXPLANATION 

CYPA 

Aerodynamic  coefficient  values 

CYPA3 

CZA 

CZA3 

Initial  condition  values 

*o 

*o 

e 

o 


Initial  condition  values 

Y 

o 

• 

Y 

o 

Z 

o 

• 

Z 

o 

MSQE  (always  0.0) 

PLOT9  data  (one  card  for  each  plot 
desired  of  i|>,  9,  or  <f>) 
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CARD 


! 


COLUMNS 

FIELD 

EXPLANATION 

1-3 

A3 

FMT(2) 

(same 

Width  of  format  for  output 
as  width  of  plot) 

1-3 

13 

LMAX: 

Width  of  plot  (1<LMAX£121) 

5-7 

13 

LO:  Value  of  initial  point 

9 

11 

LOG: 

Type  of  plot 

0  -  Linear 

1  -  Log 

2  -  Log-log 

11  A1  FORM:  Form  of  data  point  printed 

on  plot 

F  -  F  field 

E  -  E  field 

X  -  No  data  point  value  printed 
13-22  F10.3  Scale  factor 
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Case  II  Coefficient  extraction 


CARD  COLUMNS  FIELD  EXPLANATION 

1  1  II  Mode  *=  1 

2-5  Labels  for  output  (listed  at  end  of 

appendix) 

6  1-12  II  ICADJ(I)  Initial  conditions  to  be 

adjusted  l:yes,  0:no  In  order 
they  are: 


7  1-19  II  CADJ(I)  Coefficients  to  be  adjusted 

0:no,  l:yes  In  order  they  are: 

C.-,  C.-3,  C0  ,  C  C  -3,  C  C  -3,  C  -.  C  -3, 
io*  Za  Zp  mo*  ma  na’  na  npa  npa 

C  ,  C  0,  C  ,  C  -2,  C  -  C  -3  C  C  -,  C  C  -3 
mqQ  mq2  xq  xa  ya  ya  ypa  ypa  za  za 


8 

1-60 

F5.2 

QW(I)  Weight  factors  for  state  vector 
components  (same  order  as  initial  conditions) 

9 

Convergence  criteria 

1-2 

12 

MAXIT  Maximum  number  of  iterations  allowed 
before  program  is  terminated 

3-12 

F10.4 

TOL  Convergence  tolerance 

10 

1-6 

11 

IEQ(I)  0:no,  l:yes 

11 

1-12 

11 

ISV(I)  State  vector  component  -ire 
histories  to  be  read  as  data  0:no  l:yes 
(same  order  as  initial  conditions) 

12 

Integration  constants 

1-5 

F5.3 

H 

6-8 

13 

1T0 

9-13 

F5.3 

TMAX 

14-18 

F5.3 

TZERO 
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CARD  COLUMNS 


FIELD 


EXPLANATION 


19 

11 

IP  O:no  l:yes 

20 

11 

IT  O:no  l:yes 

21 

11 

IFE  0:no  l:yes 

22 

11 

NF 

constants 

1-11 

F11.8 

RO 

12-22 

F11.4 

V 

23-33 

F11.4 

AR 

34-44 

F11.4 

D 

45-55 

F11.4 

P 

Aerodynamic  constants 

1-11 

F11.8 

g 

12-22 

F11.4 

AIX 

23-33 

F11.4 

AI 

34-44 

F11.4 

AM 

Aerodynamic  coefficient  estimates 

1-10 

F10.4 

CLA 

11-20 

F10.4 

CLA3 

21-30 

F10.4 

CLP 

31-40 

F10.4 

CMA 

41-50 

F10.4 

CMA3 

51-60 

F10.4 

CNA 

61-70 

F10.4 

CNA3 

71-80 

F10.4 

CNPA 
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CARD 

COLUMNS 

FIELD 

EXPLANATION 

16 

1-10 

F10.4 

Aerodynamic  coefficient  estimates 

CNPA3 

11-2G 

F10.4 

CMQO 

21-30 

F10.4 

CMQ2 

31-40 

F10.4 

CXO 

41-50 

F10.4 

CXA2 

51-60 

F10.4 

CYA 

61-70 

F10.4 

CYA3 

71-80 

F10.4 

CYPA 

17 

1-10 

F10.4 

Aerodynamic  coefficient  estimates 

CYPA3 

11-20 

F10.4 

CZA 

21-30 

F10.4 

CZA3 

NCARDS 

1-80 

5E15.6 

Time  histories  of  state  vector  components 

18 

NCARDS 

1-10 

F10.4 

NCARDS  -NSV(NPTS/5) 

Initial  condition  estimates 

r 

11-20 

F10.4 

0 

*o 

21-30 

F10.4 

0 

31-40 

F10.4 

0 

• 

6 

41-50 

F10.4 

0 

51-60 

F10.4 

0 

• 

<b 

Yo 
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CARD  COLUMNS 


FIELD 


EXPLANATION 


61-70  F10.4  X 

o 

71-80  F10.4  X 

o 

19  Initial  condition  estimates 

NCARDS 

I- 10  F10.4  Y 

o 

II- 20  F10.4  Y 

o 

21-30  F10.4  Z 

o 

31-40  F10.4  Z 

o 

41-50  F10.4  MSQE  (always  0.0) 

20,21,22  Same  as  PL0T9  data  in  Case  I 

NCARDS 

Case  III  Flight  simulation  with  punched  output 

Same  as  Case  I  except  a  card  like  card  #11  of  Case  II  is  inserted 
between  cards  #2  and  #3  of  Case  I. 


Output  Labels  (Cards  2-5)  For  Case  II  Coefficient  Extraction 
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APPENDIX  IV 


RESULTS  OF  PROGRAM  TEST  RUNS 


ESTIMATED 

DEGREES  OF 

PARAMETER 

INITIAL 

EXTRACTED 

STANDARD 

CORRECT 

RUN 

FREEDOM 

EXTRACTED 

GUESS 

VALUE 

DEVIATION 

VALUE 

1 

1-Noise 

e 

0 

0.1800 

0.1754 

0.QG129 

0.1754* 

Linear 

• 

e 

0 

0.0 

0.0106 

0.0184 

0.0106 

C  - 
mo 

-1.750 

-2.011 

0.00541 

-2.010 

C 

mqo 

-65.000 

-60.669 

1.1834 

-60.585 

2 

1-No  Noise 

6 

o 

0.1800 

0.1745 

0.000003 

0.1745 

Linear 

• 

0 

o 

0.0 

0.0000 

0.00004 

0.0 

C  - 
ma 

-1.750 

-2.002 

0.00001 

-2.000 

C 

mqo 

-65.000 

-60.108 

0.0026 

-60.000 

3 

1-Noise 

0 

0 

0.5235 

0.5260 

0.00144 

0.5259* 

Nonlinear 

• 

0 

o 

0.0 

0.0412 

0.07249 

0.0411 

C  - 
ma 

-1.800 

-2.013 

0.00799 

-2.012 

C  -3 
ma 

-25.000 

-24.613 

0.14766 

-24.600 

C 

mqo 

-62.000 

-61.660 

0.73425 

-61.578 

Snq2 

-175.000 

-239.086 

23.947 

-238.630 

4 

1-No  Noise 

*o 

0.5235 

0.5235 

0.00003 

0.5231 

Nonlinear 

*o 

0.0 

0.0000 

0.00152 

0.0 

C  - 
ma 

-1.800 

-2.000 

0.00018 

-2.000 

C  -3 
ma 

-25.000 

-24.500 

0.00316 

-24.500 

Cmqo 

-62.000 

-60.000 

0.01551 

-60.000 

Snq2 

-175.000 

-160.004 

0.48508 

-160.000 

^Values  obtained  from  UFPLANAR 
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ESTIMATED 

DEGREES  OF 

PARAMETER 

INITIAL 

EXTRACTED 

STANDARD 

CORRECT 

FREEDOM 

EXTRACTED 

GUESS 

VALUE 

DEVIATION 

VALUE 

1-No  Noise 

0 

0 

0.5235 

0.5235 

0. 00003 

0.5235 

Nonlinear 

• 

e 

0 

0.0 

0.0000 

0.00152 

0.0 

C  - 

ma 

-1.800 

-2.000 

0.00018 

-2.000 

C  -3 
ma 

-25.000 

-24.500 

0.00316 

-24.500 

C 

mqo 

-62.000 

-60.000 

0.01548 

-60.000 

Cmq2 

-175.000 

-162.996 

0.48505 

-163.000 

2-No  Noise 

*o 

0.5235 

0.5235 

0.00001 

0.5235 

i 

0.0 

0.0000 

0.00014 

0.0 

0 

o 

0.5235 

0.5235 

0.00001 

0.5235 

• 

0 

0 

0.0 

.0000 

0.00011 

0.0 

C  - 
ma 

-1.800 

-2.000 

0.00001 

-2.000 

C 

mqo 

-62.000 

-59.998 

0.00128 

-60.000 

3-No  Noise 

0 

o 

0.5200 

0.5235 

0.00008 

0.5235 

• 

0 

o 

0.0 

0 . 0000 

0.00112 

0.0 

X 

0 

0.0 

0.0002 

0.00006 

0.0 

• 

X 

0 

550.000 

499.9993 

0 .00014 

500.000 

Z 

0 

1050.000 

999.9998 

0 .00004 

1000.000 

• 

Z 

0 

0.0 

0 .0002 

0 . 00006 

0.0 

C  - 

ma 

-1.800 

-2.000 

0 .00007 

-2.000 

C 

mqo 

-64.000 

-59.999 

0  .01661 

-60.000 

C 

xo 

0.150 

0.250 

0  .00001 

0-250 

3- No  Noise 

0 

o 

0. 1700 

0.1745 

0.00000 

0.1745 

• 

0 

0 

0.0 

0. 0000 

0.00006 

0.0 
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ESTIMATED 

DEGREES  OF 

PARAMETER 

INITIAL 

EXTRACTED 

STANDARD 

CORRECT 

RUN 

FREEDOM 

EXTRACTED 

GUESS 

VALUE 

DEVIATION 

VALUE 

X 

o 

0.0 

0.0004 

0.00025 

0.0 

• 

X 

o 

500.000 

701.1914 

0.00073 

701.19 

Z 

o 

1100.00 

999.9993 

0.00020 

1000.000 

• 

Z 

o 

0.0 

-.0010 

0.00014 

0.0 

C  - 
mo 

-11.160 

-12.399 

0.00002 

-12.400 

C 

tnqo 

-143.000 

-131.142 

0.00908 

-130.00 

C 

xo 

0.200 

0.174 

0.00051 

0.172 

Cza 

7.000 

6.553 

0.04904 

6.400 

9 

3-No  Noise 

*o 

0.4000 

0.3491 

0.00000 

0.3491 

*o 

0.0 

0.0000 

0  .00001 

0.00 

e 

o 

0 .4000 

0.3491 

0 .00001 

0 .3491 

« 

0 

0 

0.00 

0.0000 

0.00001 

0.00 

*0 

0.00 

0.0002 

0.00002 

0.00 

♦o 

0.00 

0.0011 

0 . 00009 

0.00 

0.200 

0.209 

0.00001 

0.209 

% 

-1.200 

1.292 

0.00057 

-1.325 

CmS 

-3.000 

-3.208 

0.00004 

-3.208 

c 

nq0 

-30.000 

-28.995 

0.00515 

-29.000 

c  - 

npa 

-7.000 

-7.283 

0.00219 

-7.291 

10 

2-No  Noise 

0 

o 

1.500 

1.5700 

0.00000 

1.5700 

■ 

0 

0 

0.00 

0.000 

0 . 0000 

0.00 

<t> 

ro 

0.00 

0.00 

0.00000 

0.00 

• 

<t 

y0 

0.900 

1.000 

0 . 00000 

1.000 

clp 

-1.200 

-1.000 

0.00005 

-1.000 
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ESTIMATED 

DEGREES  OF 

PARAMETER 

INITIAL 

EXTRACTED 

STANDARD 

CORRECT 

FREEDOM 

EXTRACTED 

GUESS 

VALUE 

DEVIATION 

VALUE 

C  - 
ma 

-1.800 

-2.000 

0.00000 

-2.000 

C 

mq 

-55.00 

-60.000 

0.00067 

-60.000 

-  No  Noise 

*0 

0.1700 

0.1743 

0.00004 

0.1745 

i 

0.000 

-.0001 

0.00060 

0.00 

0 

o 

0.1700 

0.1746 

0.00004 

0.1745 

• 

0 

0 

0.00 

-0.0011 

0.00057 

0.00 

♦o 

0.00 

-.0001 

0.00005 

0.00 

0.900 

1.0006 

0.00018 

1.000 

X 

o 

0.00 

0.0004 

0.00020 

0.00 

• 

X 

o 

550.00 

499.998 

0.00070 

500.00 

Y 

o 

0.00 

0.0000 

0.00015 

0.00 

• 

Y 

o 

0.00 

0.0043 

0.00028 

0.00 

Z 

o 

1100.00 

999.999 

0.00018 

1000.00 

• 

Z 

o 

0.00 

0.0066 

0.00019 

0.00 

CZp 

-1.250 

-.998 

0.00038 

-1.000 

c  - 

ma 

-1.750 

-1.990 

0.00009 

-2.000 

o 

J 

-55.000 

-59.273 

0.02063 

-60.000 

C 

a 

0 

o 

o 

CM 

• 

1 

0.249 

0.00003 

-.250 

S; 

-1.000 

-1.974 

0.00147 

-2.000 

cza 

-1.000 

-2.000 

0.00031 

-2.00C 
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APPENDIX  IV  (Continued) 


Run 

1 

2 

3 

4 

5 

6 
7 


c 


b 


c 


EXECUTION  TIME  OF  TEST  RUNS  USING 
AN  IBM  370/165  DIGITAL  COMPUTER 


Number  of 
Iterations 

4 

4 

4 

5 

5 

11 

6 
8 
6 

4 
6 
9 
2 
3 

5 

6 
7 

12 


Execution  Time 
(Seconds) 

RMSE 

18.68 

-2 

0.5484  x  10 

17.84 

0.1216  x  10-2 

30.89 

0.7027  x  10"2 

32.46 

0.1507  x  10"3 

31.56 

0.1504  x  10  3 

57.13 

0.4616  x  10  3 

51.18 

0.4259  x  10* J 

105.99 

0.1529  x  Uf  1 

/ 

0.2061  ' 

0.1578  x  \j~ 

106.91 

0.1234  >  "v 

0.1452 

0.8976  x  10* 

43.76 

0.2917 

0.7741 

215.98 

0.  5955  5  -H> 

0.2702  x 

0.1291  .. 

l 
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INITIAL  DISTRIBUTION 


ASD  (ENYS)  1 
AUL  (AUL-LSE- 70-239)  1 
COMMANDER  NWC  1 
DDC  2 
TRADOC/ADTC/DO  1 
DL  1 
DLB  1 
DLOSL  2 
USAFA/TECH  LIBRARY  1 
AFFDL  1 
USN  WEAPONS  LAB  1 
ARMY  MSL  COM/AMSMI-RDK  1 
NASA  AMES  RSCH  CENTER  1 
RSCH  LAB/ AERO  RSCH  GP  1 
USNOL  1 
AEDC/ARO,  INC.  1 
AEDC/ARO,  INC  1 
AEDC/ARO,  INC  1 
ARMY  MATERIEL  ANAL  AGCY  1 
UN IV  OF  FLORIDA  5 
DLD  1 
DLDL  2 
DLMA  15 
HQ  USAF/SAMI  1 
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